function [RMSE,MAE,R2]=evaluationIndexCalculate(y,yPredict)
%{
计算均方根误差、平均绝对误差、R方：
input:y_observed and y_predictd
output:rmse,mae,r^2
%}
n=length(y);
RMSE=sqrt(sum((y-yPredict).^2)/n);
MAE=sum(abs(y-yPredict))/n;
R2 = 1 - ( sum( (y-yPredict).^2 ) / sum( (y-mean(y)).^2 ) );

end